Intervalley Scattering and Weak Localization in Si-based Two-Dimensional Structures 
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We have measured the weak localization magnetoresistance in (OOl)-oriented Si MOS structures 
with a wide range of mobilities. For the quantitative analysis of the data, we have extended the 
theory of weak-localization corrections in the ballistic regime to the system with two equivalent 
valleys in electron spectrum. This theory describes the observed magnetoresistance and allows the 
extraction of the phase breaking time t v and the intervalley scattering time r v . The temperature 
dependences T V (T) for all studied structures are in good agreement with the theory of electron- 
electron interaction effects in two-dimensional systems. The intervalley scattering is elastic and 
rather strong: r v is typically only an order of magnitude greater than the transport time, r. It is 
found that the intervalley scattering rate is temperature- independent and the ratio r v /r decreases 
with increasing the electron density. These observations suggest that the roughness of the Si-SiC>2 
interface plays the major role in intervalley scattering. 

PACS numbers: 72.10.-d, 73.20.Fz, 73.40.Qv 
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I. INTRODUCTION 

Two-dimensional semiconductor structures with a de- 
generate ground state attract a great deal of attention 
because of the richness of low-temperature transport and 
thermodynamic phenomena that are often absent in sys- 
tems with a simple band structure. The energy spectrum 
of a two-dimensional electron system in the metal-oxide- 
semiconductor (MOS) structures grown on the (001)- 
oriented Si surface consists of six subbands (valleys). At 
low temperatures and low electron densities, only two of 
them are occupied (the other four valleys are consider- 
ably higher in energy) [lj. The low-energy valleys are 
almost equivalent: the valley splitting Ay caused by an 
asymmetry of the confining potential is typically negli- 
gible in comparison with the Fermi energy [l|. When 
the intervalley scattering is weak, the valley degeneracy 
strongly affects both electron-electron interaction and 
weak localization (WL) effects in the conductivity. In 
particular, the interaction effects in Si MOS structures 
are strongly amplified by the valley degeneracy 0; this 
accounts for the anomalous "metallic" temperature de- 
pendence of the resistivity in high-mobility Si MOSFETs 
at intermediate temperatures 0, 0] • Accordingly, the in- 
tervalley scattering plays an important role in the low- 
temperature phenomena in Si MOS structures: it de- 
termines the low-temperature cut-off of the metallic-like 
transport and could also modify the 2D metal-insulator 
transition observed in these structures at low electron 
densities [ij Q. However, to the best of our knowledge, 
there were no systematic studies of the intervalley scat- 
tering in Si MOS structures. 

The measurements of weak localization corrections to 
the conductivity of two- valley systems allow one to study 
intervalley scattering. The effect of intervalley scattering 
on WL depends on the relationship between the interval- 



ley scattering time, r v , and the time of dephasing of the 
electron wave function, t v . For weak intervalley scat- 
tering (t v >> t v ), two valleys are independent at the 
time scale t v relevant to the diffusion regime of WL. It is 
therefore expected that the magnitude of the WL correc- 
tion is doubled in comparison with its value in a system 
with strong intervalley scattering (r v << t v ); the latter 
correction is the same as in a single- valley system because 
the valleys are completely mixed at the t v time scale. 

In numerous measurements of the WL magnetoresis- 
tance in Si MOS structures @, II, H, 03 ) the experimental 
data were fitted using the Hikami-Larkin-Nagaoka (HLN) 
theory [TTj] . Interestingly, the factor-of-two enhancement 
of the WL correction was never observed, indicating that 
intervalley scattering is rather strong. In order to ex- 
tract the intervalley scattering time from the WL mag- 
netoresistance, the measurements should be extended to- 
wards higher magnetic fields. However, the HLN theory, 
which is used for fitting the WL magnetoresistance, was 
developed within the diffusive approximation (i.e small 
magnetic fields, see below). Therefore, for an adequate 
description of the effect of intervalley scattering on WL 
in Si MOS structures, a theory applicable over a wider 
range of magnetic fields should be developed. 

In this paper, we extend the ballistic (i.e. applicable to 
arbitrary classically-weak magnetic fields) theory of WL 
corrections to the case of two degenerate valleys. This en- 
ables the detailed quantitative analysis of the WL mag- 
netoresistance measured for several Si MOS structures 
with the electron mobility varying over an order of mag- 
nitude. Rather small extracted values of t v (~ 10t) indi- 
cate that the intervalley scattering in Si MOS structures 
is strong. We have found that (i) intervalley scattering 
is temperature-independent, i.e. elastic, (ii) the scatter- 
ing rate depends monotonically on the electron density, 
and (iii) there is no simple correlation between the inter- 
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valley scattering and the mobility for different samples. 
The phase relaxation time in all studied structures is well 
described by the theory of electron-electron interaction 
effects [12. 

The paper is organized as follows. In section [II] we de- 
scribe the samples, measurement techniques and provide 
examples of experimental data. In section UTIl to compare 
our data with previously reported results, we fit the WL 
magnetoresistance with the HLN theory. The results of 
this analysis indicate that the WL correction is approx- 
imately a factor of 2 smaller than one could expect for 
a system with two independent valleys. The theory of 
the WL magnetoconductance for systems with degener- 
ate multi- valley spectrum is presented in section IIVI and 
the data are analyzed using this theory in section [Vj In 
section PVTl we discuss the results of the data analysis, in 
particular, the intervalley and phase breaking times. 



II. EXPERIMENT 

Below we present data for three representative Si MOS- 
FET samples: (i) high-mobility (fi » 2 m 2 /Vs at O.IK) 
sample Si6-14 [J] which demonstrates a strongly pro- 
nounced "metallic-like" dependence cr(T) and a metal- 
insulator transition with decreasing electron density n, 
(ii) sample Si39 with an intermediate mobility \i w 
0.45 m 2 /Vs and (iii) low-mobility sample Si40 with 
fi « 0.18 m 2 /Vs at T < 4.2 K. The transport times 
for these samples within the studied range of n were 
t pa 2 ps, 0.6 ps, and 0.2 ps, respectively. The resistance 
was measured using the standard four-terminal AC tech- 
nique by a resistance bridge LR700 (for sample Si6-14) 
and by SR830 lock-in amplifier in combination with a 
differential preamplifier (for samples Si39 and Si40). A 
sufficiently small measuring current (1-3 nA for Si6-14 
and 50-100 nA for Si39 and Si40) was chosen to avoid 
overheating of electrons. The electron density, controlled 
by the gate voltage, was found from the period of the 
Shubnikov-de Haas oscillations and the dependence of the 
Hall resistance on magnetic field. Both results were con- 
sistent with each other within 2% accuracy in the studied 
range of densities; this uncertainty of n is insignificant for 
the further analysis. To ensure that only the lowest size 
quantization subband is filled [l[ , we performed measure- 
ments at n < 4 x 10 12 cm~ 2 . 

The WL magnetoresistance was measured at T = 
0.05 — 0.6 K for high-mobility sample Si6-14, and at 
T = 1.3 - 4.2 K for samples Si39 and Si40. At these tem- 
peratures, the phase breaking time exceeds the transport 
time by one to two orders of magnitude. The magnetic 
field aligned perpendicular to the plane of Si MOS struc- 
tures was varied from -1 to +lkG (Si6-14) and from -3 
to +3kG (Si39 and Si40). When sample Si6-14 was mea- 
sured at T < IK, an additional in-plane field ~200G 
was applied to quench the superconductivity in the cur- 
rent/voltage contact pads and the gate electrode made of 
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FIG. 1: Examples of the magnetoresistance p(B) data for Si6- 
14, n = 1 X 10 12 cm _2 (a), Si39, n = 3.6 x 10 12 cm -2 (b), and 
Si40, n = 3.5 x 10 12 cm" 2 (c). The data within the hatched 
regions have been used to extract the t v value. Upper axes 
show the magnetic field in units of Btr = <E"o /27rZ 2 , lower axes 
show the field in Gauss. The data points were chopped for 
clarity. 



thin aluminum film. For reliable extraction of the phase 
relaxation time t v from the WL magnetoresistance, we 
have chosen a small field step size: 1 G for Si6-14 and 3 G 
for Si39 and Si40. The examples of magnetoresistance 
p{B) data for Si6-14, Si39 and Si40 at a fixed density and 
various temperatures are shown in Figs. QJl, [TJd and[]J;, 
respectively. Hereafter throughout the paper we will use 
magnetoconductance (MC) Act = p(B)^ 1 — p(B = 0) _1 . 



III. FITTING THE DATA WITH THE 
HIKAMI-LARKIN-NAGAOKA THEORY 



It is a common practice to extract the phase breaking 
time from the WL magnetoconductance using the HLN 
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FIG. 2: a) Example of the MC data (points) for sample Si40, 
T = 1A5K and n =3.34xl0 12 cm" 2 . Solid line is calculated 
using Eq.(l) with two fitting parameters a — 1 and t v /t — 
133. The dashed curve is an attempt to fit the same data 
over the same range of b with a fixed prefactor a = 2 and 
t^/t =61. b) Dependences of the fitting parameters a (solid 
squares) and t>/t (open squares) on the magnetic field range 
0-6 which was used for fitting. 
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tem for various values of r^/r and magnetic fields b. It 
was found that both approaches agree with each other 
within a limited range of fields b < 0.15 provided that 
T ip/ T > 30 and the conductivity is much greater than 
e 2 /2ir 2 h. Thus, within these limits, Eq. |T]) can be used 
for extraction of t v from the WL magnetoresistance in 
a single-valley system, and the adjustable parameter a 
is aproximately equal to 1. For a system with two val- 
leys and weak intervalley scattering, the prefactor a is 
expected to be two times larger, because each valley con- 
tributes the term Acthln with a ~ 1 to Ac. 

Figure [3]a shows a fit of our typical WL MC curve with 
Eq. {T]). The fitting performed over the magnetic field 
range b = - 0.2 gives t^/t = 133 and, contrary to the 
expectation for a two-valley system without intervalley 
scattering, a — 1. Changing the magnetic field range, 
where the data are fitted, causes only minor variations of 
these parameters (see Fig. [2]b). An attempt to analyze 
the MC curve using Eq.([T]) with a fixed prefactor a = 2 
results in a much worse fit (dashed line in Fig. [2]a) . 

Several reasons for the reduction of a in a single valley 
system have been considered in Ref. [HI], including (i) 
Maki-Thompson correction, (ii) Density-of-states correc- 
tion, (iii) higher order corrections in {k^l) -1 , where Uf 
is the Fermi wave vector, and (iv) low t v /t ratio. The 
corrections (i) and (ii) were shown to be small [14| . To 
ensure that the higher order corrections are also small, 
we have studied the WL MC only for large conductances 
~ 100 x e 2 /2n 2 Ti. When the ratio t v /t decreases, the 
HLN theory becomes inadequate, and the fitting proce- 
dure results in an artificially reduced prefactor. Corre- 
spondingly, we performed measurements at such low tem- 
peratures that the inequality t v /t > 30 was satisfied. We 
conclude therefore that the aforementioned reasons can- 
not account for a low value of the prefactor a ~ 1 in the 
studied multi-valley structures. It should be noted that 
tt « 1 fa Si MOS structures was obtained in numerous 
previous experiments 0, d, H, [l(| • We show below that 
the prefactor reduction can be well described by the the- 
ory which explicitly takes the intervalley scattering into 
account. 



Here if> is the digamma-function, e is the electron charge, 
Ti is the Planck constant, b = B/Bt r is the dimensionless 
magnetic field, B tr — Q> /2irl 2 , <E> = nh/e, and I is the 
transport mean free path [l5j]. The prefactor a and the 
dimensionless ratio t v /t are treated as fitting parame- 
ters. Note that with an increase of the magnetic field, 
the crossover from the diffusive regime (b << 1) to the 
ballistic regime (b ~ 1) is expected in the WL correc- 
tions. Equation |T]) with prefactor a = 1 is the exact 
result for a single-valley system in the diffusive regime, 
i.e. at 7^ 3> r and for sufficiently small magnetic fields 
b <C 1 11]. On the other hand, the experimental data are 
often obtained beyond these limits and, therefore, should 
be described by more general ballistic theories [3, [13] • 

In Ref. [HI], the HLN theory was numerically com- 
pared with the ballistic theory for a single-valley sys- 



IV. THEORY 

A consistent theory of weak localization is developed 
in the framework of the diagram technique. The weak- 
localization corrections to the conductivity arise in the 
first order in the parameter (kpl) , where the mean 
free path I is governed by the scattering time r, which is 
controlled by both intra- valley (r v ) and inter- valley (Ti) 
scattering processes: 



i/t = i/t v + i/n 



(2) 



The weak-localization correction to the conductivity in 
the magnetic field has the form 



Atr(S) = A(t (q) + Aa w 



(3) 
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where the terms Aa^ and Aa^ correspond to the 
standard dia gra ms, which have been considered in detail 
in Refs. [13, 13, H(| • We neglect both valley and spin- 
orbit splitting in Si MOS structures. [2l[ Then, for the 
short-range scattering potential [22| , one obtains 



Aff(a) = £E / J 2 Mp')cT a a \p',P)d P d P > 

a/3 



(4) 



af3-yS 

x [G«(p, p')G fl (p", p) + G A (p, p')G A (p", p)]dpdp'dp" • 

Here a, /3,7, (5 = 1,2 are the valley indices; Ja;(p, p') is 
the ir-component of the current vertex, which is defined 
as 



E W Z J p (p, p")c ( ^{ P ", p 1 ) dp" , 
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where F(p,p') = G A (p, p')G R (p, p>). 

To calculate the weak-localization correction to the 
conductivity in multi-valley structures, we assume that 
the impurity potential is the same for particles in differ- 
ent valleys and, thus, the electron scattering in the valleys 
is strongly correlated. In particular, in the (OOI)-oriented 
Si-based structures the nonzero correlators are 
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J(p,p') 



ie^f—^[G A (p,p>) + G*(p,p>)}, (5) 

TO* |p — p'| 



G A (p, p') and G fl (p, p') are the advanced and retarded 
Green functions, 
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^,ky(pW s ,kSp') 



E F - E s ± ih/(2r) ± ih/(2T v 



(6) 



^Ps,k (p) is the wave function of an electron subject to 
an external magnetic field B, which is perpendicular to 
the 2D plane, s and k y are the quantum numbers (s is 
the Landau level number and k y is the in-plane wave vec- 
tor), E s = huj c (s + 1/2) is the energy of the s th Landau 
level, LOc = eB/m* is the cyclotron frequency, to* is the 



effective electron mass; C,!^ (p, p') and C. 



(3)a/3, 



- 7(5> (P,P')are 

the Cooperons which depend on four valley indices. The 



parameters W^j? are determined by intervalley and in- 
travalley correlators of the scattering potential and are 
defined by 



Other correlators vanish due to averaging over the im- 
purity positions because the lowest conduction-band val- 
leys in silicon are located in the A-points of the Brillouin 
zone and, therefore, the Bloch functions contain oscilla- 
tory factors. We note that the intravalley and interval- 
ley scattering times are determined by these correlators: 
l/ Ti = m*WH/h 3 , 1/tv = m*W&/h 3 . 

Using the standard procedure (see Refs. [l|J[23[), one 
can expand P(p, p') and the Cooperons in the series 
of eigenf unctions of a particle with the double charge 
in a magnetic field and derive equations for the weak- 
localization corrections Aa^ and Ad' 6 ' . Calculations 
show that the corrections have the form 



2 7 oo 
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Aff(b) = ^5* E ^ + C ^ Qn/ 2 ■ (H) 

N=0 



(V Q fe 1 , ( 3fc 2 T / 7 fc3,5fc 4 )Mmp = W"f <5fc 1+ fe 3 ,fe 2+ fe 4 , (7) 

where V a k 1 ,/3k 2 is t ne matrix element of scattering be- 
tween electron states (/3, k%) an d (a, k\) in zero magnetic 
field, kj (j = 1 ... 4) are wave vectors in the 2D plane, 
-^Vimp is the two-dimensional density of impurities, and 
the angle brackets stand for the averaging over impurity 
spatial distribution. 

The Cooperons C^ a/3 (p, p') and C^ Q/3 (p, p') repre- 
sent the sums of internal parts of the fan diagrams start- 
ing with two and three lines, respectively, [13, EH l20j |. 
They can be found from the following equations: 



C$ a/3 (p,p')= > 
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2(1 - t /t v ) 3 P n , P N (1 - 2t/t v ) 3 P j 



/v 



1 - (1 - t/t v )P n 1-P n 1 - (1 - 2t/t v )P n 

(12) 

where P/v and Qn are coefficients which are given by 
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and Ijn and Ljy- are the Laguerre polynomials. 
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FIG. 3: a) WL magnetoconductance calculated for a two- 
valley system using Eqs. I3I10I11I (solid lines) and for a single- 
valley system [l7l ] (dotted line). For solid curves from bottom 
to top, t v /t = 10, 25, 50, 100, 1000, 10000, respectively. Two 
upper curves are indistinguishable by eye. For all the curves 
t v /t = 100. The hatched region was used for fitting with 
Eq. |T}. b) Dependences of the fitting parameters, a and 
t v /t, on t v /t. The data points at r v /r = 1 correspond to a 
single- valley system. 



Equations ©, pTJj) . and (fTTjl describe the WL magne- 
toconductance over the whole range of classically weak 
magnetic fields uj c t = \iB < 1. In the limit of vanishing 
intervalley scattering (l/r v = 0), Eqs. (fTTTjl and (fTTj) are 
reduced to the conventional expressions for the WL cor- 
rections to the conductivity of a single- valley system [ijj 
and, in particular, to the HLN formula in the diffu- 
sion regime. The only difference is a prefactor of 2, which 
accounts for the valley degeneracy. 

To illustrate the effect of intervalley scattering on the 
magnetoconductance, we calculated the Ao~(b) depen- 
dence using Eqs. ©, 1[ID]). and (JTTJ) for a fixed tJ/t = 100 
and various values of t v /t. The results are shown in 
Fig. [3]a by solid lines. For comparison, we also calcu- 
lated the MC using a similar theory jl7| developed for a 
single- valley system (dotted line). We then fitted these 
dependences over the range b < 0.15 using the HLN the- 
ory [Eq. (1)] with two fitting parameters, the prefactor a 
and t v /t. In other words, we fitted the theoretical curve 
the same way as the experimental data have been fitted 
above in Sec 

Figure [3]b shows the resultant fitting parameters; for 



completeness, we also depicted a and r^/r for a single- 
valley system at t v /t = 1. The main results of the fit 
are as follows: (i) the extracted phase breaking time t v 
coincides with its true value within a few percent (this 
uncertainty is insignificant for further analysis), and (ii) 
the observed prefactor increases from w 1 to « 2 as t v 
increases and becomes greater than t v . Therefore, the 
approximate equality as; 1 is simply a consequence of a 
large ratio t v /t v 3> 1. 



V. FITTING THE DATA WITH THE 
BALLISTIC THEORY 

It is intuitively clear and will be discussed in more 
detail below, that the MC in low fields b <C 1 is predomi- 
nantly determined by t v . In principle, one could use the 
"ballistic" theory for fitting the MC data in the whole 
range of magnetic fields and thus for determining both 
t v and t v from a single fit. However the series Eqs. (fTUj) 
and (|11[) converge very slowly in small b region. There- 
fore, to determine t v , it is more practical to use in low 
fields Eq. (Q]) , that was shown (see Fig. [3]) to provide the 
correct t v value. 

Consequently, we have used the following procedure of 
extracting r v from the WL magnetoconductance. Firstly, 
we analyzed the MR data in sufficiently weak magnetic 
fields and at low temperatures. In this regime (the 
hatched regions in Fig. [J), the dephasing occurs at a time 
scale much greater than r v , and we can apply Eq. ([1]) for 
extracting t v ; the second adjustable parameter, prefactor 
a, appears to be close to 1. At the next stage, we substi- 
tute t v into the "ballistic" formulae [Eqs. ©, ([TO]). (fPT|) ] 
and calculate the MC curves in a wide range of fields 
(b < 1) for various r v . Figure 0] illustrates this procedure 
using as an example the same MC data as in Fig. [5] We 
calculate Aer (b) using the summation technique similar 
to that described in Ref. [24[ for single- valley systems. 

Figure 0] shows that the experimental MC (circles) is 
smaller than the MC for a system with two unmixed de- 
generate valleys (curve 1) and larger than MC for a single 
valley system (curve 5) . This observation again indicates 
that the MC in the studied Si MOS structures is affected 
by valley mixing. Curves 2, 3, and 4 in Fig. 2| correspond 
to t v /t = 15, 12, and 9, respectively. Note that in the 
magnetic field range b < 0.15, these three curves, the 
experimental data, the HLN formula, and the ballistic 
result for a single valley system are almost indistinguish- 
able from each other (see the inset to Fig. Q|. Therefore, 
t v cannot be reliably found from the MC in low fields. 

Figure [4] shows that the discrepancy between the curve 
5 for a single- valley system and the curves 2,3,4 f° r tw° 
mixed valleys (r v /r =15, 12, and 9) grows as b increases. 
This observation has a transparent physical explanation: 
with increasing b, the typical size of electron trajectories 
which contribute to the WL correction diminishes, and 
the valley mixing over the time of travel along these tra- 
jectories becomes small when b > &q/Dt v . As a result, 



6 
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FIG. 4: Comparison between the WL magnetoresistance for 
sample Si40, n = 33.4 x 10 u cm~ 2 , T = 1.45 K and the 
ballistic theory. Different MC curves are calculated using 
Eqs. (l3)), (|10[) , (|ll[) : 1 - for two unmixed valleys (r v = oo); 
2 - t v /t =15; 3 - t v /t = 12; 4 - t v /t = 9; 5 - MC for 
a single- valley system (equal to curve 1 divided by 2). In- 
set blows up the data in the range b < 0.15. Curve 6 is the 
HLN theory (Eq. JTJ|) with a prefactor a = 1 (see Fig. [2]a). 
t v /t = 133 for all calculated curves. 



the WL magneto-conductance in strong magnetic fields 
approaches the theoretical prediction for a two- valley sys- 
tem with no intervalley scattering. 

We also note that all calculated curves deviate from 
the experimental data. As Fig. 0] shows, curve 4 calcu- 
lated for t v /t = 9 at 6 > 0.4 is approximately paral- 
lel to but lower than the experimental data in magnetic 
fields 6 > 0.4. On the other hand, curve 2 calculated for 
t v /t = 15 almost coincides with the data in low magnetic 
fields b < 0.4, though deviates substantially from them 
in higher fields. The minimal mean-square deviation of 
the calculated curve from the data is realized for r v /r 
=12 (curve 3). 

Thus, the value of r v depends on the magnetic field 
interval (61,62) within which the MC data is fitted. 
The (t v /t) values, obtained from fitting the difference 
Act(&i) — Aer(6 2 ) as a function of (61,62), decrease as 6 = 
(61 + 6 2 )/2 increases (see Fig. [5]). This monotonic depen- 
dence has been reproduced for all samples and tempera- 
tures. We believe that this apparent r v (6) dependence is 
an artifact of the fitting procedure. In all above calcula- 
tions we assumed t v to be field- independent. However, 
t v should depend on a perpendicular magnetic field [12] • 
To the best of our knowledge, there are neither experi- 
mental nor theoretical systematic studies of this depen- 
dence beyond the diffusive limit. Ignoring this depen- 
dence in our fitting could lead to the observed monotonic 
variation in r v with 6. 
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FIG. 5: Intervalley scattering time determined from fitting 
the difference, Acr(bi) — A<x(&2), for the same magnetoresis- 
tance curve as in Fig. [4] The fitting ranges (61, 62) are shown 
by the horizontal bars. t v /t = 133. 

A question arises therefore what range of magnetic 
fields should be chosen for t v extraction? To answer this 
question we have analyzed errors of our method; the re- 
sulting root-mean-square sum of all errors is shown by the 
error bars in Fig. O The error analysis is presented in de- 
tail in the Appendix (Sec. IIX[) , where it is shown that nei- 
ther small fields (6 < 0.1) nor large fields (6 ~ 1) should 
be used for r v extraction. In weak fields, the MC is insen- 
sitive to r v , whereas in strong fields one approaches the 
limits of applicability of the theory developed in Sec. IIVI 

Therefore, we conclude that an intermediate range of 
magnetic fields is most suitable for extracting r v . For the 
further analysis, we choose the range 6 = 0.2 - 0.4. We 
have verified that our conclusions on the temperature and 
density dependences of r v are not affected if this range is 
changed. 



VI. RESULTS AND DISCUSSION 

A. Phase breaking time 

As we have already mentioned, at the first stage of 
the analysis we estimated the phase breaking time r^. 
Comparison of the t v (T) dependences with the theory 
of interaction effects [12[ is shown in Figs. [5] a-c. The 
uncertainty in the values of a and t v (shown as the 
error bars in Fig. [6]) reflects mainly the uncertainty in 
a(b) in the weak fields 6 < 0.01. The magnitude of 
the phase breaking time and its temperature dependence 
are in good agreement with the theory for all samples 
within the studied ranges of electron density [Si6-14: 
n = (0.28 - 1.5) x 10 12 cm" 2 , Si39: n = (2 - 2.5) x 10 12 
cm -2 , Si40: n = (3 - 4) x 10 12 cm" 2 ]. Note that no 
adjustable parameters are involved in this comparison, 
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FIG. 6: Temperature dependence of the extracted t v value 
in units of r (left axes) and in picoseconds (right axes): a) 
Si6-14, n = 9.98xlO n cm~ 2 , b) Si39, n = 29.4x 10 n cm- 2 , 
c) Si40, n = 33.4xlO n cm- 2 . Solid lines show the t v {T) 
dependence predicted by the theory of interaction corrections 
[12 ] with 15 triplet terms, dashed line - with 3 triplet terms. 
The insets show the corresponding temperature dependences 
of the prefactor a. 



5% (dashed line in Fig. [6]). The condition Ti/t v < ksT 
is violated at temperatures lower than 0.3 K for sample 
Si6-14. Still, the t v {T) data agree better with theory 
when all 15 rather than 3 triplet terms are taken into 
account. The observed quantitative agreement of the ex- 
perimental values of t v with the theory suggests that 
t v is weakly affected by intervallcy scattering near the 
crossover ksTT v /h ~ 1. 



B. Prefactor a 

By fitting the weak-field MC data with the HLN the- 
ory, we obtained the prefactor a that is close to 1 for 
all samples (see the insets to Fig. [6J; this suggests that 
the valleys are intermixed on the t v time scale. The 
decrease of a from 0.9 to 0.6 with increasing tempera- 
ture, obtained for sample Si39 (see the inset to Fig. 6b), 
we believe, is an artifact, because relatively small values 
t v /t ~ 30, observed for this sample at high tempera- 
tures, make Eq. (1) inadequate. The complete theory 
described in Sec. IIVI explains that small value of a is 
a consequence of a fast phase relaxation. For the same 
reason, there is a larger scattering in the values of a(T) 
for samples Si6-14 and Si40 at the highest temperatures 
(Fig. |BJ), where r v is small. 

It is worth noting that the smallness of prefactor a has 
been attributed to the intervalley scattering in Ref. 
However, the MC data in this experiment were fitted 
with the theory [23|, which does not take into account 
the non-backscattering correction Eq. (fTTj) . Our esti- 
mates show that for the parameters of samples studied 
in Ref. [8] (t v /t = 20, r v /r = 4, and b = 0.5-7), the 
non-backscattering correction contributes about 50% to 
the extracted value of r v . 



Intervalley scattering time: independence of 
temperature 



the Fermi-liquid parameter Fq and the effective mass 
m* were obtained in independent measurements [25| . 

The theoretical curves (solid lines in Fig. \§§ are calcu- 
lated following Ref. [l2| for 15 triplet channels ||, which 
implies small valley splitting (fc^T > A v ) and relatively 
weak intervalley scattering (H/t v < fcgT). As found from 
the analysis of the low-temperature transport and mag- 
netotransport data [26|], the condition fceT > A v was sat- 
isfied for samples Si6-14 and Si39 over the major part of 
the studied temperature range. Whether or not this con- 
dition is fulfilled for low mobility samples Si39 and Si40 is 
actually not important, because the measurements were 
performed at such high densities that the amplitude of 
the triplet term in the interaction corrections to t v was 
small in comparison with the singlet term: changing the 
number of triplet terms from 15 (two-valley case) to 3 
(single-valley case) caused variation of t v by less than 



Following the procedure described in Sec. El we have 
extracted r v by fitting the WL MC data with "ballistic" 
Eqs. ((5|).([ICJ | .((TT] 1 . Figure [7] shows that the values of t v 
are temperature-independent within the accuracy of our 
measurements. This observation suggests that the inter- 
valley scattering is elastic, i.e. governed by static disor- 
der. Similar conclusion can be also drawn from the fact 
that a remains close to 1 , while the extracted t v much ex- 
ceeds r v and grows without saturation as T decreases (see 
Fig. [6]). Indeed, were the intervalley scattering inelastic, 
one would have observed a prefactor a ~ 2 because the 
dephasing would occur in two valleys independently and 
the intervalley scattering would be just an additional de- 
phasing mechanism. In the latter case, a cut-off of the 
dephasing time at the level t v — r v is also expected. The 
two observations, the absence of the cut-off and a ~ 1 
support the self-consistency of our analysis. 
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FIG. 7: Temperature dependence of tv in units of r (left axes) 
and in picoseconds (right axes): a) Si6-14 n=9.98xl0 11 cm~ 2 , 
b) Si39 n=29.4xl0 n cm- 2 , c) Si40 n=33.4x 10 n cm~ 2 . Solid 
horizontal lines show the average r v . 



The intervalley transitions are not expected to be in- 
elastic for the following reason. The intervalley scattering 
requires a large momentum transfer comparable to the 
vector of reciprocal lattice 2ir ja ~ 10 8 cm _1 (a is the in- 
teratomic distance). At liquid helium temperatures only 
static disorder can cause these transitions, as the mo- 
menta of electrons kp ~ 10 6 cm -1 for the studied range 
of densities and phonons k p h ~ fc^T '/ (hv s ) (here v s is the 
sound velocity) arc much smaller than 2ix j a. Static disor- 
der can lead only to elastic scattering because it changes 
momentum of scattered electrons but does not change 
their energy. 



D. Intervalley scattering time: density and sample 
dependence 

Figure [8] shows the density dependence of r v values 
averaged over the temperature. For all three samples, the 
relative rate of the intervalley transitions (with respect to 
the momentum relaxation rate) increases with density. 
This points to the dominant role of the Si-SiC>2 interface 
in the intervalley transitions. The electron wave function 

in Si MOS structure is positioned mostly in the bulk 
silicon and exponentially decays in SiC>2 [1]. When the 
gate voltage (and, hence, the density n) is increased, the 
electrons are "pushed" towards the Si-SiC>2 interface, and 
the amplitude of the wavefunction at the interface, \I/o ; 
increases. The probability of the interface scattering is 




Density (10 11 cm" 2 ' 



FIG. 8: Density dependence of the intervalley scattering time 
(averaged over temperature) for samples Si6-14 (a), Si39 (b) 
and Si40 (c). 



proportional to l^ol 2 an d increases with n [l[; this is in 
line with the behavior shown in Fig. [5] 

In the experiments we used samples with the mobilities 
which vary over a decade. We find no correlation between 
t v and the mobility for different samples. This suggests 
that the intervalley scattering is determined by a sample- 
specific interface disorder, namely the surface roughness 
at the atomic length scale, which might be different for 
the samples fabricated on different wafers. In contrast to 
the intervalley scattering, the mobility is governed mostly 
by impurities in the bulk and by the interface roughness 
at a large length scale, ~ 2ir/k~F. 

The measured values of t v for all samples are within 
the interval (3-12)r, which indicates that the valley index 
remains a good quantum number at the time scale ~ r. 



VII. SUMMARY 

To summarize, we have studied the weak localization 
magnetoconductance in Si MOS structures over wide 
ranges of the electron densities, mobilities, and tempera- 
tures. In order to quantitatively analyze the experimen- 
tal data, we have developed the theory of weak localiza- 
tion for two-dimensional multivalley systems, which is 
valid in both the diffusion and ballistic regimes. The the- 
ory, which explicitly takes the intervalley scattering into 
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account, allowed us to conduct the first detailed study of 
the intervalley scattering in the Si MOS structures. It 
was found that: 

1. Intervalley scattering in Si MOS structures is an 
elastic and temperature-independent process. 

2. The ratio t v /t monotonically increases as the elec- 
tron density decreases. This observation suggests 
that the intervalley scattering is governed by the 
disorder at the Si-Si02 interface. 

3. There is no simple correlation between the inter- 
valley scattering rate and the sample mobility (or 
the momentum relaxation rate); this points to a 
sample-specific rather than universal mechanism of 
the intervalley scattering. 

4. The smallness of the prefactor a ~ 1, that is ob- 
tained from fitting the experimental WL data with 
the HLN formula, is a consequence of a fast in- 
tervalley relaxation rate which exceeds the phase 
relaxation rate. 

5. The temperature dependence of the phase relax- 
ation time in Si MOS structures is in quantitative 
agreement with the theory of electron-electron in- 
teraction effects in disordered two-dimensional sys- 
tems [Hj]. 

We note that the approach similar to that developed 
in this paper can be used for studies of intervalley re- 
laxation in other multi-valley two-dimensional electron 
systems, such as AlAs-AlGaAs heterostructures 12711 . Si 
MOX structures (28|, and Si-SiGe quantum wells 29]. 
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IX. APPENDIX: ANALYSIS OF POSSIBLE 
ERRORS IN tv 

We present here an analysis of errors in the fitting pro- 
cedure which determine the size of error bars in Fig. [S] 
As discussed in Section [V] r v was found from the follow- 
ing equation: 



Act(6i)exp - Act(6 2 )exp = (14) 

= Aa(bi,T v /r, t v /t) T h - Act(6 2 , t v /t, t v /t)th 

Here subscript EXP denotes experimental data, subscript 
TH denotes calculation using Eqs. (|3")l . (fTU)) . (fTTj) . Conse- 
quently, the uncertainty in t v is determined by (i) uncer- 
tainty in b, (ii) uncertainty in t v and (iii) uncertainty in 
the conductivity. To estimate each contribution to the 
error, we varied the corresponding parameter (6, t v or 
Act) within its uncertainty and determined the variation 
in r v by solving Eq. (jT4"|) . 

The uncertainty 8b in b = 2Birl 2 /§o value is deter- 
mined by the uncertainty in the mean free path I [l5|. 
The latter is about 2-3 % due to the uncertainties in 
electron density n and Drudc conductivity. However, 8b 
affects t v rather weakly for the following reason: MC 
in the studied magnetic field range behaves approxi- 
mately as ln(fe), therefore Act(&i) — Act(62) ~ hi(&i/&2) = 
\n{B x /B 2 ). 

The error related to the uncertainty in t v is essential 
in low magnetic fields where magnetoconductance is sen- 
sitive to t v . Correspondingly, the error bars in low fields 
b < 0.15 in Fig. [5] are determined predominantly by the 
uncertainty in t v . 

Another source of errors is related to the precision of 
the absolute value of WL MC ( "calibration error" ) . In- 
deed, the accuracy of our measurements of the absolute 
magnetoresistance value is ~ 0.5%. Higher order cor- 
rections, Maki-Thompson and DOS corrections (lij can 
modify MC by approximately 2-3% (as shown in Ref. [l4| , 
£(Act)/Act « 2e 2 pu/2TT 2 h « 0.025). In order to estimate 
this error we artificially changed our experimental data 
by 3% and studied the corresponding change in r v . The 
error appears to grow in small magnetic field where mag- 
netoconductance is weakly sensitive to t v . Therefore, 
small fields should not be used for the extraction ofr v . In 
large magnetic fields (6 ~ 1) MC becomes again weakly 
sensitive to r v , and the latter error grows as b increases, 
as shown by the error bars in Fig. The calibration 
error is minimal in intermediate magnetic fields, where 
MC is most sensitive to r v . 

Our attempts to analyze the WL MC data in strong 
fields b > 1 using Eqs. (|3"J).([TU ]) .(fTT ]) resulted in a large 
uncertainty of the fitting parameter t v /t (large scatter- 
ing of extracted t v /t for various electron densities and 
temperatures). In large magnetic fields, there are sev- 
eral other error mechanisms which are difficult to take 
into account. For example, at b ~ 1, t v differs from 
its small- field value Moreover, in Ref. [H the MC 
for b > 1 was shown to behave in a non-universal man- 
ner: it strongly depends on details of scattering poten- 
tial, whereas our theory assumes an uncorrelated short- 
range disorder. Some other mechanisms of magnetocon- 
ductance (such as classical memory effects, interaction 
corrections, Maki-Thompson corrections etc.) may also 
become essential in large fields where the shape of WL 
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MC curve flattens. Therefore we believe that the in- extraction of the intervalley scattering rate, 
termediate field range b — 0.2 — 0.4 is optimal for the 
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